General unit-disk representation for periodic multilayers 
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We suggest a geometrical framework to discuss periodic layered structures in the unit disk. The 
band gaps appear when the point representing the system approaches the unit circle. We show that 
the trace of the matrix describing the basic period allows for a classification in three families of 
orbits with quite different properties. The laws of convergence of the iterates to the unit circle can 
be then considered as universal features of the reflection. 



Photonic band gap structures [lj] can be dealt with by 
computing the band structure (using e.g. Bloch theory) 
or from the point of view of scattering 2\. The essential 
difference is that a scattering experiment always involves 
a finite structure, while Bloch waves imply an endlessly 
repetition of the basic period. 

In the context of electromagnetic optics, photonic crys- 
tals (that is, one-dimensional periodic layered structures) 
have attracted recently a lot of attention because the 
striking property of acting as omnidirectional reflectors: 
they can reflect light at any polarization, any incidence 
angle, and over a wide range of wavelengths 0, El 0> B • 

The appearance of strong reflection (stop bands) de- 
pends on the properties of the basic period. Different 
theoretical approaches, involving equivalent medium the- 
ories, group velocity, and other quantities ja, have 
been developed to account for the detailed structure of 
these stop bands and their edges. Each one of these mod- 
els emphasizes some aspects of the problem but, at the 
same time, has some drawbacks. 

In the present Letter we introduce a geometrical set- 
ting that allows for a deeper understanding of periodic 
systems. Our treatment is quite general and only as- 
sumes linearity: it applies to any physical system whose 
transfer matrix belongs to the group SU(1, 1). The key 
point for our purposes is the fact that the multilayer 
transfer function induces a bilinear transformation in the 
unit disk [Tol ITl| . Since perfect mirrors are represented 
precisely by the unit circle, the route to a stop band can 
be understood as the convergence of the point represent- 
ing the action of the system to the unit circle. 

We start by examining the basic period of the system, 
which consists of a stack of plane-parallel layers sand- 
wiched between two semi-infinite ambient (a) and sub- 
strate (s) media that we shall assume to be identical, 
since this is the common experimental case. Hereafter 
all the media are supposed to be lossless, homogeneous, 
and isotropic. 

A monochromatic linearly polarized plane wave falls 
from the ambient making an angle Oq with the normal 

to the first interface and with an amplitude E a + \ We 
consider as well another plane wave of the same frequency 
and polarization, and with amplitude Eg , incident from 
the substrate at the same angle 6*o- The output fields in 
the ambient and the substrate will be denoted E a ' and 
E^ + \ respectively. 

The field amplitudes at each side of the multilayer are 



related by the linear relation 
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Here the complex numbers R as and T as are, respec- 
tively, the overall reflection and transmission coefficients 
for a wave incident from the ambient. Because |i? as | 2 + 
\T as \ 2 = 1, we have det M as = +1 and then the set of loss- 
less multilayer matrices reduces to the group SU(1, 1). 

We are often interested in the transformation proper- 
ties of field quotients rather than the fields themselves. 
Therefore, we introduce the complex numbers 
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for both ambient and substrate. Equation 0) defines 
then a transformation on the complex plane C, mapping 
the point z s into the point z a , according to 



Za = 3>[M as 



P* + a*z s 
a + (3z s 



(4) 



which is a bilinear (or Mobius) transformation. This ac- 
tion can be seen as a function z a = f(z s ) that will be 
called the multilayer transfer function. One can check 
that the unit disk, the external region and the unit circle 
remain invariant under the multilayer action |10|. Note 
that \z a \ = \z s \ — 1 for a perfect mirror with light in- 
cident from both the ambient and the substrate, so the 
transformation J3J relates then points on the unit cir- 
cle. When no light strikes from the substrate z s = and 
\z a \ = 1, so a mirror maps the origin into a point on the 
unit circle. 

In what follows, the idea of fixed points of the trans- 
formation will prove to be essential. These points can be 
defined as the field configurations such that z a = z s = zt 
in Eq. (@J; i.e., Zf = &[M as ,Zf], whose solutions are 
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± { -2t Im(a) ± V[Tr(M as )P - 4} . (5) 
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FIG. 1: Plot of typical orbits in the unit disk for: (a) canonical 
transfer matrices as given in Eq. and (b) arbitrary transfer 
matrices. 



The trace of M as provides then a suitable tool for the 
classification of multilayers [l2j . 

When [Tr(M a ,,)] 2 < 4 the multilayer action is elliptic 
and it has only one fixed point inside the unit disk, while 
the other lies outside. When [Tr(M as )] 2 > 4 the action is 
hyperbolic and it has two fixed points on the unit circle. 
Finally, when [Tr(M as )] 2 = 4 the multilayer action is 
parabolic and it has only one (double) fixed point on the 
unit circle. 

To proceed further let us note that by taking the con- 
jugate of M as with any matrix C £ SU(1, 1), that is 
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we obtain another matrix of the same type, since 
Tr(M as ) = Tr(M as ). Conversely, if two multilayer ma- 
trices have the same trace, a matrix C satisfying Eq. © 
can be always found. 

The fixed points of M os are then the image by C of the 
fixed points of M as . In consequence, given any matrix 
M as , it can always be reduced to a unique M as with one 
of the following canonical forms: 



A(x) 



exp(i(p/2) 

exp(— iip/2) 

cosh(x/2) i sinh(x/2) 
—i sinh(x/2) cosh(x/2) 

l-iry/2 77/2 
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that have as fixed points the origin (elliptic), +i and 
—i (hyperbolic) and +i (parabolic), and whose physical 
significance has been studied before |l5| . 

The concept of orbit is especially appropriate for ob- 
taining a picture of these actions. Given a point z, its 
orbit is the set of points z' obtained from z by the action 
of all the elements of the group. In Fig. l.a we have plot- 
ted typical orbits for each one of the canonical forms. For 
matrices K(ip) the orbits are circumferences centered at 
the origin. For A(x), they are arcs of circumference going 
from the point +i to the point —i through z. Finally, for 
N(?7) the orbits are circumferences passing through the 
point +i and joining the points z and — z* . 




Elliptic 



Hyperbolic 



Parabolic 



FIG. 2: Plot of the successive iterates (N — 1, . . . , 5) for an 
elliptic, hyperbolic, and parabolic action starting from the 
origin as the initial point. The physical system is described 
in the text. Only hyperbolic and parabolic actions tend to 
the unit circle. 



For an arbitrary matrix M as the orbits can be obtained 
by transforming the orbits described before with the ap- 
propriate matrix C . In Fig. l.b we have plotted typi- 
cal examples of such orbits for elliptic, hyperbolic, and 
parabolic actions. We stress that once the fixed points of 
the matrix M as are known, one can ensure that z a will 
lie in the orbit associated to z s . 

Assume now that we have a finite periodic structure 
obtained by repeating N times the basic period repre- 
sented by M as . The overall transfer matrix for this sys- 
tem is M^. In the unit-disk picture, the transformed field 
by the A^-period structure is represented by the point 



Z N = $[M oa ,2jV-l] = *[M£,J8b] 



(8) 



where zq is the initial point z s . The idea that maps 
iterates could be applied to this problem was recognized 
before, though with a somewhat different approach [Tfi[ 

m 

Henceforth, we shall take zq = 0, which is not a seri- 
ous restriction and corresponds to the case in which no 
light incides from the substrate [E& ' = 0] , as it happens 
usually. Note also that all the points z n lie in the or- 
bit associated to the initial point zq by the single period, 
which is determined by its fixed points: the character 
of these fixed points determine thus the behavior of the 
periodic structure. 

To illustrate how this geometrical approach works in 
practice, we take the single period as a Fabry-Perot-like 
system formed by two identical plates separated by a 
spacer of phase thickness 62 ■ This is a symmetric system 
for which R as and T as can be easily computed By 
varying 82 we can choose to work in the elliptic, the hy- 
perbolic, or the parabolic case. In Fig. 2 we have plotted 
the sequence of successive iterates obtained numerically 
for these three regimes. 

In the elliptic case, it is clear that the points z jv revolve 
in the orbit centered at the fixed point and the system 
never reaches the unit circle. 

On the contrary, for the hyperbolic and parabolic cases 
the iterates converge to one of the fixed points on the 
unit circle, although with different laws, which corre- 
spond to the band stop and band edges of the system, 
respectively Q. 

To gain further insights into these behaviors, we com- 
pute explicitly the ATth iterate. This can be easily done 
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for the canonical forms in Eq. (JTJ and then, conjugating 
as in © we obtain, after some lengthy calculations, that 
for a hyperbolic action one has 



iN 



ZN 



l-Z»(Zf+/Zf-) 



(9) 



where £ = (a + f3zf-)/(a + /3z/+) is a complex number 
satisfying |£| < 1. Analogously, for the parabolic case we 
have 



ZN 



N(3z 2 



f 



Nfizf - 1' 



(10) 



where Zf is the (double) fixed point. It is quite obvious 
that in both cases zn converges to one of the fixed points 
on the unit circle, so \zn\ — > 1 when N increases, a typical 
behavior of perfect mirror. In the mathematical litera- 
ture this limit point is referred to as the Denjoy- Wolff 
point of the map 0] . 

To characterize the convergence of Zn we note that, 
because zq = z s = 0, this initial point is transformed by 
the single period into z a — R as . Therefore, zn represents 
the reflection coefficient of the overall periodic structure 

i2da , which is obviously different from (R as ) N . One can 



then compute that for the hyperbolic case |l£ 



ZN 



[sinh(x)/ sinh(iVx)p 



(11) 



that approaches the unit circle exponentially with N, 
as one could expect from a band stop, while for the 
parabolic action 



ZN 



| 2 + (1/A0 2 ' 



(12) 



that goes to unity with a typical behavior 0(N~ 2 ). This 
is universal in the physics of reflection, as put forward in 
a different framework by Yeh [2(j and Lekner [2l| . 

To conclude, we expect that the geometrical scenario 
presented here could provide an appropriate tool for an- 
alyzing and classifying the performance of periodic mul- 
tilayers in an elegant and concise way that, additionally, 
is wider enough to accommodate other periodic systems 
appearing in physics. 

We acknowledge Jose F. Carinena and Jose M. Mon- 
tesinos for illuminating discussions. 
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